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Abstract 

We  discuss  some  issues  of  displaying  two-dimensional  surfaces  in  four-dimensional 
(4D)  space,  including  the  behavior  of  surface  normals  under  projection,  the  silhouette 
points  due  to  the  projection,  and  methods  for  object  orientation  and  projection  center 
specification.  We  have  implemented  an  interactive  4D  display  system  on  z-buffer  based 
graphics  workstations.  Preliminary  experiments  show  that  such  a  4D  display  system 
can  give  valuable  insights  into  high-dimensional  geometry.  We  present  some  pictures 
from  the  examples  using  high-dimensional  geometry,  offset  curve  tracing  and  collision 
detection  problems,  and  explain  some  of  the  insights  they  convey. 

1  Introduction 

The  geometry  of  high-dimensional  space  has  shown  to  be  quite  useful  in  the  area  of  CAGD 
and  solid  modeling.  Applications  include  describing  the  motion  of  3D  objects,  modeling 
solids  with  nonuniform  material  properties,  and  formulating  constraints  for  offset  surfaces 
and  Voronoi  surfaces  [3,  8, 15].  For  inspection  and  understanding  of  the  properties  related  to 
geometry,  pictures  are  very  effective  to  provide  intuitions.  Unfortunately  it  is  very  hard,  if 
not  impossible,  for  us  to  visualize  objects  in  high-dimensional  space.  Therefore,  visualization 
of  high-dimensional  space  by  means  of  computer  graphics,  especially  interactively,  is  a 
research  topic  that  attracts  growing  attention  [10,  15]. 

‘Work  supported  in  part  by  NSF  grant  CCR  86-19817  and  DMC  88-07550  and  by  the  office  of  Naval 
Research  under  contract  N00014-90-J-1599. 
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Visualization  of  4D  space  is  a  good  starting  point  because  not  only  is  it  relatively 
easy,  but  also  many  problems  fit  naturally  into  a  4D  formulation.  For  example,  describing 
3D  objects  in  motion,  embedding  3D  projective  space  into  affine  space,  and  analyzing 
plane  curves  with  complex  roots  and/or  coefficients.  Since  4D  space  has  many  properties 
unfamiliar  to  us,  its  visualization  is  full  of  problems  to  be  explored. 

We  briefly  review  some  basic  ideas  useful  in  4D  visualization.  A  solid  object  in  4D  space 
is  of  dimension  four.  Its  boundary  is  composed  of  one  or  more  hypersurfaces  of  dimension 
three.  In  nondegenerate  cases,  the  intersection  of  two  hypersurfaces  is  a  surface  of  dimension 
two,  and  the  intersection  of  three  hypersurfaces  is  a  curve  of  dimension  one.  In  contrast  to 
3D  space,  two  surfaces  in  4D  space  generally  intersect  in  a  point  instead  of  a  curve.  Since 
human  beings  have  no  sense  of  4D  space,  we  have  to  map  the  4D  solid  from  4D  space  into  3D 
space.  One  method  is  to  intersect  the  4D  object  with  a  hypersurface  (perhaps  a  hyperplane 
normal  to  one  coordinate  axis)  and  get  an  image  in  3D  space.  In  order  to  perceive  the  whole 
4D  object,  a  series  of  images  with  different  positions  of  the  intersecting  hypersurface  would 
be  needed.  Another  method  is  to  project  the  4D  object  into  3D  space,  orthographically  or 
perspectively.  Again,  a  series  of  images  with  different  directions  of  projection  is  needed. 
Both  methods  must  be  supplemented  with  one  more  intersection  or  projection  step  if  we 
want  to  render  the  3D  image  on  a  conventional  2D  device  such  as  a  piece  of  paper  or  a 
computer  screen.  The  intermediate  step  of  mapping  objects  into  3D  can  be  so  implemented, 
or  can  be  combined  with  the  mapping  to  2D  into  a  single  procedure. 

Eckhart  [4]  proposes  a  method  to  project  an  object  into  several  planes  that  are  orthog¬ 
onal  to  different  pairs  of  coordinate  axes.  The  2D  images  so  obtained  can  be  put  together 
in  a  systematic  way  in  analogy  to  the  principal  views  in  traditional  engineering  drawings. 
However,  since  this  method  only  displays  curves  that  are  three  dimensions  lower  than  a  4D 
object,  and  since  the  viewer  must  gather  information  from  different  pictures,  it  may  be  very 
hard  to  interpret  such  pictures. 

This  paper  concentrates  on  the  interactive  display  of  two  dimensional  surfaces  in  4D 
space,  thereafter  referred  to  as  2-surfares.  It  is  a  worthwhile  job  for  several  reasons: 

1.  Some  problems  are  naturally  formulated  with  2-surfaces  in  4D  space. 

2.  2-surface  display  is  an  important  subtask  in  visualizing  4D  objects. 

3.  Interactive  speeds  are  possible  because  of  the  available  hardware  on  graphics  devices. 
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For  example,  consider  the  collision  detection  problem  of  two  solids  moving  in  3D  space. 
Their  moving  boundaries  are  hypersurfaces  in  4D  space.  The  intersection  of  the  two  hy¬ 
persurfaces  is  generally  a  2-surface  in  4D  space.  Also,  for  any  curve  problem  with  three 
equations  in  four  variables,  we  can  think  of  two  of  the  equations  as  the  definition  of  a  2- 
surface,  and  consider  the  third  equation  as  the  constraint  for  a  point  of  the  2-surface  to 
be  on  the  curve.  Displaying  the  curve  and  the  2-surface  simultaneously  is  usually  easier  to 
understand  than  displaying  the  curve  alone. 

2-surface  display  is  an  important  subtask  in  the  visualization  of  4D  solids.  A  4D  solid  can 
be  displayed  by  its  boundaries  which  are  hypersurfaces.  The  projection  of  a  hypersurface 
into  3D  space  is  a  3D  volume.  Basically  we  have  three  choices  in  displaying  the  3D  volume: 
volume  rendering,  surface  rendering,  and  curve  rendering  (wire  mesh).  In  volume  rendering, 
each  voxel1  corresponds  to  several  points  on  the  hypersurface.  It  can  be  reduced  into  one 
point  on  the  hypersurface  if  we  use  a  hidden  hypersurface  removal  technique,  similar  to  the 
hidden  surface  removal  technique  used  in  3D  rendering.  To  obtain  a  nontrivial  picture,  we 
need  a  shading  model  in  4D  space  [15].  Although  this  is  an  interesting  topic,  we  may  want 
to  eliminate  the  effect  of  4D  light  while  trying  to  understand  the  structure  of  hypersurfaces 
as  a  set  of  points.  On  the  other  hand,  rendering  a  4D  solid  by  curves  shows  a  wire  mesh, 
which  involves  no  shading  model  at  all.  It  is  frequently  used  for  displaying  surfaces  in  3D 
space.  However,  it  seems  hard  to  interpret  such  pictures  when  they  are  used  for  displaying 
hypersurfaces  of  a  4D  solid  [13].  Rendering  a  4D  object  via  surface  rendering  shows  a  group 
of  surfaces  placed  within  a  hypersurface,  and  serves  as  a  bridge  between  volume  rendering 
and  curve  rendering.  Surface  rendering  involves  only  a  3D  shading  model  [7].  Here,  the  role 
of  the  2-surfaces  for  displaying  a  hypersurface  is  analogous  to  that  of  curves  for  displaying 
a  surface  in  3D  space.  Often  we  can  gain  valuable  insights  from  this  analogy. 

Interactive  display  is  crucial  for  4D  visualization.  After  projection  from  4D  to  3D,  and 
finally  to  2D,  a  significant  amount  of  information  is  lost.  To  compensate,  the  viewer  should 
be  able  to  see  a  real  time  animation  of  the  object  in  translation  and  rotation,  controlled 
by,  say,  a  mouse.  Currently  most  graphics  workstations  are  designed  for  3D  rendering,  and 
use  efficient  techniques  of  surface  display.  Investigating  of  the  relationship  between  surfaces 
before  and  after  projection  may  lead  to  efficient  algorithms  for  4D  visualization. 

We  are  experimenting  with  an  interactive  4D  visualization  system  on  conventional  z- 

lA  voxel  is  a  three-dimensional  pixel 
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buffer  based  graphics  workstations.  The  input  2-suifaces  are  first  polygonalized  by  algebraic 
or  space  division  methods  [1].  The  polygons  in  4D  space  are  projected  into  3D  space  and 
then  fed  into  the  3D  graphics  engine.  The  method  of  “polygonalization  before  projection” 
is  more  desirable  than  the  method  of  “projection  before  polygonalization”  for  interactive 
display.  Usually,  polygonalization  requires  more  computation.  Therefore,  polygonalizing 
the  2-surface  as  a  preprocessing  step  means  that  a  better  response  can  be  obtained  when 
changing  the  projection  parameters  repeatedly.  Moreover,  polygonalization  in  4D  space 
can  better  account  for  the  intrinsic  geometric  properties  of  the  2-surface.  Some  of  these 
properties  are  distorted  by  the  projection  to  3D  space.  We  have  carefully  designed  a  user 
interface  which  offers  position  and  orientation  control  of  objects,  lights,  and  projection 
centers. 

The  paper  is  organized  as  follows.  Section  2  deals  with  intrinsic  geometric  issues  such 
as  surface  normals  and  silhouettes  under  projections.  Section  3  discusses  one  aspect  of 
the  user  interface,  namely,  how  to  extend  the  Euler  angles  into  4D  space  and  use  these 
angles  uniformly  for  object  orientation  and  projection  direction  specification.  In  section  4, 
several  examples  are  discussed  and  illustrated  with  pictures.  Section  5,  finally,  draws  some 
conclusions  from  this  work. 

2  Surface  Normals  And  Silhouettes  Under  Projection 

We  map  a  2-surface  in  4D  first  to  a  surface  in  3D  space.  Then  the  surface  is  rendered  on 
a  2D  device  by  the  standard  3D  methods  for  shading  and  illumination.  To  determine  the 
intensity  of  the  light  reflected  by  a  surface  in  3D  space,  it  is  necessary  to  find  the  non.ial 
of  the  surface  at  each  point.  However,  the  projections  of  the  normals  of  a  2-surface  in 
4D  do  not  necessarily  coincide  with  the  normal  of  the  projected  2-surface.  0  ie  way  to 
find  the  normal  of  the  projected  2-surface  is  to  calculate  it  from  the  equation  representing 
the  projected  2-surface.  Another  way  is  to  calculate  the  normal  directly  from  the  normal 
space  of  the  2-surface  before  projection.  The  latter  seems  to  be  mor'  efficient  because  the 
derivation  of  the  equation  of  the  projected  2-surface  could  be  expensive  [9]. 

In  this  section  we  discuss  in  detail  how  to  calculate  the  nor.ial  of  the  projected  2-surface 
from  the  normals  of  the  2-surface  before  projection.  Given  a  2-surface  and  a  projection, 
there  are  certain  points  at  which  the  normal  space  Hoes  not  determine  the  normal  of  the 


4 


projected  2-surface.  We  will  show  that  they  are  exactly  those  points  whose  tangent  space 
is  projected  into  a  space  of  lower  dimension,  and  that  they  are  silhouettes.  We  also  discuss 
how  to  recognize  silhouettes. 

A  2-surface  5  in  4D  space  ('ll4)  can  be  represented  either  in  implicit  form: 

f\(x,y,z,w )  =  0 
f2(x,y,z,w)  =  0 

or  in  parametric  form: 

x  =  fi(s,t)  z  = /3(M) 
y=f2(s,t)  w  =  f4(s,  t) 

At  a  nonsingular  point  p  on  S,  there  is  a  unique  normal  space  and  a  unique  tangent 
space.  The  two  spaces  are  orthogonal  to  each  other.  Since  in  H4  both  normal  space  and 
tangent  space  of  a  2-surface  have  dimension  two,  both  are  determined  either  by  two  linearly 
independent  normal  vectors  ni  and  n2,  or  by  two  linearly  independent  tangent  vectors  ti 
and  t2.  For  a  2-surface  in  implicit  form,  it  is  easy  to  find  the  two  normal  vectors  at  a 
nonsingular  point  p: 

ni  =  V/i(P) 
n2  =  V/a(p) 

On  the  other  hand,  for  a  2-surface  in  parametric  form,  it  is  easy  to  find  the  two  tangent 
vectors  at  a  point  p  if  the  corresponding  parameters  s  and  t  are  given: 

.  _  ,  dfi  df2  dfz  df4  T 

1  -  1  da  ’  ds  ’  ds  ’  ds  ’ 

_  ,  dfi  df2  dfi  dfj  T 

2  _  1  dt'  dt'  dt'  dt  ’ 

A  simple  method  can  be  used  to  find  two  tangent  vectors  ti  and  t2  from  two  normal  vectors 
ni  and  n2  and  vice  versa.  Let  i,  j,  k,  1  be  the  base  vectors  of  TI4,  and  a,  b,  c  be  three 
vectors  where  a  =  (ax,av,az,aw)T  and  so  on.  We  define  operation  0  as: 

i  j  k  1 

ax  ay  az  aw 

®(a,b,c)  = 

bx  by  bz  bw 

C  J-  Cy  cz  CW 
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Obviously  ®(a,b,c)  is  orthogonal  to  the  subspace  span(a,b,c)  if  a,b,c  are  linearly  in¬ 
dependent.  Given  ni  and  112,  choose  any  two  vectors  a  and  b  such  that  nj,n2,a,b  are 
linearly  independent.  A  base  of  the  tangent  space  is  then 


ti  =  ®(ni,n2,a) 
t2  =  ®(ni,n2,b) 


In  most  cases  a,  b  can  be  chosen  from  the  base  vectors,  and  doing  so  reduces  the  needed 
arithmetic  operations.  Sometimes  b  can  be  replaced  by  ti  to  make  ti  and  t2  orthogonal 
to  each  other. 

A  projective  mapping  v  from  TZ4  to  TZ3  can  be  written  as 


»(p)  =  v(Pr,Py,Pz,Pw)  =  (  r-_P*  •»  T~~~  ,  7~  f 

1  rp%u  1  ^Pw  1  rpw 

where  we  assume  the  center  of  projection  is  at  I/r  =  (0,0,0,l/r)r.  It  is  an  orthographic 
projection  for  r  =  0  and  a  perspective  projection  for  r  >  0.  The  projection  center  can  be 
put  into  an  arbitrary  position  if  translation  and  rotation  in  TZ4  are  done  before  applying  v. 

Given  a  2-surface  S  and  its  normal  plane  N  at  a  point  p,  by  projection  we  get  their 
images  in  TZ3  denoted  by  v(S),  v(N),  and  r(p)  respectively.  In  general  v(N)  is  of  dimension 
two  and  cannot  be  used  as  the  normal  of  v(S)  at  point  v(p).  Fortunately  tangency  is  a 
projectively  invariant  property  (see,  e.g.  [16]).  Hence,  a  natural  method  to  calculate  the 
normal  of  the  projected  surface  is  the  following: 


Step  1:  From  the  two  normals,  find  two  tangent  vectors  tj  and  t2  as  explained  above. 


Step  2:  Project  the  three  points  p,  p  ■+■  tj,  p  +  t2  into  TZ3  with  the  mapping  v. 


Step  3:  Use  the  cross  product  to  find  the  normal  n 

n  =  [t>(p  +  ti)  -  v(p)]  x  [v(p  +  t2)  -  v(p)] 


To  derive  a  more  efficient  method,  we  need  to  know  the  relationship  of  the  2-surface 
normals  before  and  after  projection.  Given  a  vector  a  =  (ax,ay,az,zw)T  in  TZ4,  we  define 
its  natural  projection  into  TZ3  as  JrxYz(a)  =  (or, ay, az)T.  Given  a  vector  a  =  ( ax,ay,az)T 
in  TZ3,  we  define  its  natural  extrusion  into  TZ4  as  £w(a)  =  (ar,ay,a**0)r. 

Lemma  1  Suppose  that  ti,t2  (ni,n2)  are  two  linearly  independent  tangent  (normal) 
vectors  at  a  nonsingular  point  p  on  a  2-surface  5,  that  v  is  a  projection  with  center  1/ r , 
and  that  r  =  1/r  -  p  is  the  ray  from  the  point  p  to  the  projection  center.  Then, 
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(a)  if  aj  =  £w(^(p  +  tj)  —  t>(p)),  then  aj,t;,r  are  linearly  dependent  for  i  =  1,2. 

(b)  if  n  =  ®(ti,t2,r),  then  JTxYz(n)  is  the  normal  vector  of  v(S)  at  point  v(p). 

(c)  if  n  -  ani  +  /?n2  satisfying  n  •  r  =  0,  then  xxYz(n)  is  the  normal  vector  of  v(S)  at 

point  v(p). 


Proof:  (a)  Since  there  is  no  difference  for  i  =  1,2  we  drop  the  subscripts  of  t;  and  aj 
temporarily.  The  vector  a  can  be  written  as: 


a 


Then 


PT  +  tx 

Px 

l-r{pw  +  tw) 

l-rpu, 

Pv+‘y 

Py 

l_r(Pui+*u;) 

1  — rpu, 

P,  +  t, 

Pt 

l-rfpw  +  iu,) 

1  ~rpw 

V  o  ) 


aja  +  a2t  +  03(l/r  -  p)  =  0 


where 


«i  =  (1  -  rp,,,  -  rtw)(l  -  rpw) 

a2  =  — (1  —  rpw) 

O3  =  Ttw 

Clearly  a2  is  nonzero  because  otherwise  p  is  mapped  to  infinity.  Moreover,  by  choosing  the 
length  of  t  appropriately,  oj  is  also  nonzero.  For  orthographic  projection,  we  can  simply 
set  r  =  1  and  03  =  tw,  and  the  proposition  is  still  true. 

(b)  Let  a;  (i  =  1,2)  be  defined  as  in  (a).  Then, 

*XYz(n)  •  [v(p  +  tj)  -  v(p)]  =  ®(t!,t2,r)  •  aj 

=  det([ai,ti,t2,r]r) 

The  determinant  is  zero  according  to  (a). 

(c)  n  =  ani  +  /3n2  guarantees  that  n  •  ti  =  0  and  n  •  t2  =  0.  Together  with  n  •  r  =  0 

we  know  that  n  is  in  the  same  direction  as  ®(ti,t2,r).  By  (b)  the  conclusion  follows.  □ 

When  the  normal  vectors  of  a  surface  are  directly  available,  the  method  described  in 
part  (c)  of  the  lemma  is  quite  efficient.  This  would  be  the  case  for  definitions  of  the  form 
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f(x>  J/t^w)  =  0  g(x,y,z,w).  When  the  tangent  vectors  are  directly  available,  the  method 

described  in  part  (b)  can  be  used.  This  would  apply  for  parametric  surface  definitions. 

Lemma  1  cannot  be  applied  to  calculate  the  normal  in  all  cases.  If  r  lies  in  span(tj,  t2), 
the  calculated  n  will  be  a  zero  vector.  This  case  cannot  be  avoided  with  a  different  choice 
of  tj,t2  or  ni,ri2,  and  is  a  property  of  the  projection. 

We  define  silhouette  points  on  a  surface  5  with  respect  to  a  projection  v  as  those  non¬ 
singular  points  whose  tangent  space  reduces  dimension  under  the  projection  v. 

Lemma  2  Suppose  that  ti,t2  are  two  linearly  independent  tangent  vectors  at  point  p  on 
a  surface,  that  v  is  a  projection  with  center  1/r,  and  that  r  =  1/r  -  p  is  the  ray  from  p  to 
the  projection  center.  Then  p  is  a  silhouette  point  with  respect  to  v  if  and  only  if  ti,t2,r 
are  linearly  dependent. 


Proof:  First  we  assume  that  ti,t2,r  are  linearly  dependent,  and  so  r  can  be  expressed 
as  linear  combination  of  tj,t2-  Let  aj  (i  =  1,2)  be  the  natural  extrusion  of  the  vector 
w(p  +  tj)  -  v(p).  By  Lemma  1  (a)  aj,  tj,r  are  linearly  dependent.  Then, 


aj  =  <*iti  +  02^2 

*2  =  /Ml  +  Ait  2 


Consider  the  fourth  coordinates,  ai„,  and  a2tlp,  of  the  above  equations.  Then  aju,  =  a 2W  =  0 
by  definition.  But  and  1 2u,  cannot  be  both  zero  because  rw  is  not  zero.  So 


<*1  «2 
A  A 


=  o 


Hence  aj,a2  are  in  the  same  direction  and  so  are  their  natural  projections.  The  direction 
will  not  change  if  we  choose  any  other  pair  of  vectors  within  the  tangent  plane,  which  means 
the  tangent  plane  becomes  a  line  under  projection.  So  p  is  a  silhouette  point. 

Now  assume  ti,t2,r  are  linearly  independent.  Again  by  Lemma  1  (a)  we  have 


*1  =  cMi  +  a2r 

&2  =  /M2  +  Ar 

Since  both  ai  and  A  are  nonzero,  ai,a2  cannot  be  linearly  dependent.  The  tangent  plane 
does  not  reduce  dimension  under  projection.  So  p  is  not  a  silhouette  point.  □ 
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The  picture  we  actually  see  is  two-dimensional.  So  we  are  more  concerned  with  the 
silhouette  points  with  respect  to  a  projection  from  TZ4  to  TZ2.  Define  v  as: 


«(p)  =  v{px,py,pz,pw)  =  (• 


Pr 


Py 


- ; - 

1  -nPw-  r2pz  ’  1  -  Tipw  -  r2pz 

where  1/rj  and  k/r2  are  the  two  projection  centers.  Note  that  the  projection  order  is 
irrelevant:  We  get  the  same  result  mapping  first  from  TZ4  to  the  XYZ-hyperplane  and  then 
to  the  XY-plane,  as  we  obtain  by  mapping  first  from  TZ4  to  the  XYW- hyperplane  and  then 
to  the  XY-plane. 


Lemma  3  Suppose  that  ti,t2  (ni,n2)  are  two  linearly  independent  tangent  (normal) 
vectors  at  a  nonsingular  point  p  on  a  2-surface  S,  that  v  is  a  projection  from  TZ4  to  TZ 2  with 
the  centers  1/rj  and  k/r2,  and  that  ri  =  1/rx  -  p,  r2  =  k/r2  -  p  are  the  two  rays  from  the 
point  p  to  the  two  projection  centers.  Then, 

(a)  if  a;  =  £zw(f(p  +  tj)  ~  u(p))t  then  ab*iirl»r2  are  linearly  dependent  for  i  -  1,2. 

(b)  p  is  a  silhouette  point  with  respect  to  v  if  and  only  if  ti,t2,ri,r2  are  linearly  depen¬ 

dent. 


(c)  let  mi,  m2  be  two  linearly  independent  vectors  in  the  plane  orthogonal  to  the  plane 
span(ri,r2),  then  p  is  a  silhouette  point  if  and  only  if  ni,n2,mi,rri2  are  linearly 
dependent. 

Proof:  (a)  Since  there  is  no  difference  for  i  =  1,2  we  drop  the  subscripts  of  tj  and  a; 
temporarily. 

^  _ , _ Pat**  - - Em -  ^ 

l-n(Pui  +  <u>)-ra(p,+t»)  l-riPw-rjp, 

_  Pv  +  ty  _  _ Pfl _ 

a  _  l-n(pti,  +  t<i,)-r3(p,+t,)  l-riPuz-rap, 

0 

\  0 

Then, 

aia  +  a2 1  +  03ri  +  o4r2  =  0 


where 


o i  =  [1  -  r\{pw  +  <„,)  -  r2(pz  +  t*)](l  -  ripw  -  T2pz) 

o2  =  -(1  -  riPw  -  r2Pz) 
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<*3  =  ritw  +  rxr2(tzpw  -  twp2) 

^4  =  ^2 tz  ^"1^*2 (.^iPw  ~~  t\vpz) 

Clearly,  c*2  is  nonzero  because  otherwise  p  is  projected  to  infinity.  Moreover,  by  choosing 
the  length  of  t  appropriately,  ai  is  also  nonzero. 

(b)  First  assume  that  ti,t2,ri,r2  are  linearly  dependent.  Since  span(tj,t2)  is  of  di¬ 
mension  two,  we  can  find  another  vector  b  such  that  rj,r2  each  is  a  linear  combination  of 
ti,t2,b.  By  (a)  we  have 


»1  =  aiti  +  £*2*2  +  03b 

«2  =  +  ^2^2  +  03^ 

That  means  ai,a2  lie  in  the  intersection  of  span(ti,t2,b)  and  span(i,j).  This  intersection 
must  be  of  dimension  lower  than  two  because  otherwise  span(ti,  t2,  b)  is  equal  to  span(i,  j,  c) 
for  some  vector  c,  which  causes  rj,r2  to  be  parallel  within  ZW-plane,  and  so  p  is  mapped 
into  infinity.  Consequently  p  is  a  silhouette  point. 

Now  assume  that  t1,t2,ri,r2  are  linearly  independent.  Again  by  (a)  we  have 

ai  =  axti  +  «2ri  +  a3r2 
»2  =  /M2  +  $2  Fl  +  fa,*2 

Since  both  a\  and  /S;  are  nonzero,  ai,a2  must  be  linearly  independent,  and  so  are  their 
natural  projections  into  V? .  Therefore  p  is  not  a  silhouette  point. 

(c)  It  is  sufficient  to  prove  that  ti,t2,ri,r2  are  linearly  independent  if  and  only  if 
ni,  it2, mi, m2  are  linearly  independent.  Assuming  that  tx,t2,ri,r2  are  linearly  indepen¬ 
dent,  they  form  a  base  of  TZ4.  Then  both  ni  and  112  are  linear  combinations  of  ri,r2, 
and  both  mi  and  m2  are  linear  combinations  of  ti,t2-  Therefore  ni,n2,mi,m2  must  be 
linearly  independent.  The  converse  direction  is  symmetric.  □ 

We  will  explain  the  application  of  the  above  lemma  by  the  examples  in  Section  4. 
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3  Orientation  Specification 

3.1  Object  orientation  specification 

Rigid  body  rotation  in  4D  space  can  be  expressed  by  an  orthonormal  4x4  matrix  A  called 
the  direction  cosine  matrix: 

P  =  4Pl 

In  the  equation,  p  =  ( x,y,z,w)T  is  a  vector  expressed  in  the  world  coordinate  system  and 
Pi  =  (xi,yuZi,wi)T  is  the  same  vector  expressed  in  the  body-fixed  coordinate  sysiem. 
Among  the  sixteen  elements  in  the  direction  cosine  matrix  only  six  are  independent.  It  will 
be  convenient  to  specify  the  orientation  of  an  object  by  six  independent  parameters.  One 
way  is  to  write  the  direction  cosine  matrix  as  a  product  of  six  basic  rotation  matrices.  A 
basic  rotation  is  a  one-parameter  rotation  within  a  plane  spanned  by  two  coordinate  axes. 
For  example, 

'  1  0  0  0 

0  cos  0  —  sin  0  0 

#yz(0)  = 

0  sin  0  cos  0  0 

^00  0  1 

is  a  rotation  within  the  plane  span(j,k)  whose  normal  plane  is  span(i,l).  The  six  basic 
rotations  should  be  chosen  systematically  such  that  the  geometric  relationship  is  easy  to 
explain  and  easy  to  remember.  We  have  extended  into  4D  the  Euler  angles  commonly  used 
in  3D  rigid-body  kinematics  [17]. 

3D  Euler  angles  specify  the  orientation  of  objects  as  three  successive  rotations,  first 
about  the  z-axis,  then  about  the  x-axis,  and  finally  about  the  z-axis  again.  Note  that  these 
axes  are  body-fixed.  If  the  3D  Euler  angles  are  01,02*03,  the  relationship  between  the  world 
coordinates  p  and  the  body-fixed  coordinates  p3  is  given  by: 

p  =  #Xy(01  )7?Yz[02)^Xy(03)P3 

We  conceptualize  Euler  angles  as  two  separate  rotation  phases:  In  the  first  phase,  spec¬ 
ified  by  f?x y(0i )#Yz(02),  the  body-fixed  z-axis  is  oriented  in  H3  aid  put  into  its  final 
position.  In  the  second  phase,  a  single  rotation  in  7 Z2  orthogonal  to  the  already  oriented 
body-fixed  z-axis  brings  the  entire  object  into  its  final  orientation.  Thus  the  second  phase 
is  expressed  as  Rxy{0^). 
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To  clearly  distinguish  the  two  phases,  we  add  a  number  subscript  to  the  rotations, 
indicating  the  subspace  that  is  rotated.  Thus  the  coordinate  relationship  equation  by  3D 
Euler  angles  can  be  rewritten  as: 

P  =  #XY3(0l)#YZ3(02)#XY2(03)P3 

Note  that  conceptually  Rx Y3  is  a  rotation  in  7l3  and  Rx Y2  is  a  rotation  in  7 Z2,  but  their 
matrix  forms  are  the  same  if  both  are  written  as  an  n  x  n  matrix,  for  n  >  2. 

We  extend  Euler  angles  into  7Z 4  as  follows: 

1.  Orient  the  body-fixed  w-axis  of  the  object  by  three  rotations  in  1Z4: 

RXY4(#1  )R Y24(02  )# ZW4(03) 

2.  Orient  the  body-fixed  z-axis  of  the  object  in  the  3D  subspace  orthogonal  to  the  already 
oriented  tn-axis,  by  two  rotations  in  R3:  RxYziO^R-yizi^s)- 

3.  Orient  the  remaining  body-fixed  x-  and  y- axes  in  the  2D  subspace  orthogonal  to  the 
already  oriented  z-  and  tu-axes,  by  a  single  rotation  in  R2:  RxY2(9e)- 

It  is  clear  that  with  this  conceptualization  Euler  angle  specification  can  be  naturally 
extended  to  any  dimension.  The  Euler  angle  specification  in  dimension  n  +  1  is  done 
inductively  by: 

1.  Orienting  the  new  xn+1-axis  using  n  angles. 

2.  Orienting  the  n-dimensional  subspace  orthogonal  to  the  x,.+i-axis  with  the  n-dimensional 
Euler  angle  specification. 

Note  that  the  two  dimensional  Euler  angle  is  simply  a  single  rotation  within  a  plane. 

In  4D  space,  the  relationship  between  the  world  coordinates  p  and  the  body-fixed  coor¬ 
dinates  pe  is: 

P  =  £xY4(0l)#YZ4(02)tfzW4(03)-ftxY3(04)#YZ3(05)#XY2(06)P6 

The  rotations  caused  by  each  phase  are  shown  in  Figure  1.  In  the  following  equations 
Pi  =  (xi,  yi,  Zi,  Wi)T  is  the  same  vector  expressed  in  the  i-th  coordinates  system,  that  coin¬ 
cides  with  the  body-fixed  coordinate  system  at  the  end  of  the  i-th  rotation.  Particularly, 

P3 <  P5 1 P6  are  the  coordinates  at  the  end  of  the  three  rotation  phases,  respectively. 
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—  RxY4(0i)Ry?4(02)R2W4(03) 


Cl  -S1C2  Sl52C3  --Sia2S3 
•Si  CiC2  -C1S2C3  Ci3233 


0  S2 

0  0 


-C2S3 


=  -SxY3(04)-RyZ3(05) 


=  RxY2^6) 


C4  ->S4C5  S4S5  0 

S4  C4C5  -C4S5  0 

0  s5  c5  0 

0  0  0  1 

Cg  -  Sg  0  0  \  f  x 
eg  0  0  y 

U  0  1  0  2 

0  0  0  1  xl 


The  direction  cosine  matrix  in  7£4  can  also  be  considered  as  the  world  coordinates  of 
the  four  base  vectors  i,j,k,l  of  the  body-fixed  coordinate  system  at  the  end  of  the  rotation: 

*X  jx  kx  lx 
ly  jy  ky  ly 

A  = 

iz  jz  kz  lz 

\  jw  kw  lw  J 

It  is  easy  to  verify  that  the  corresponding  Euler  angles  can  be  found  as  follows: 

Step  1:  Calculate  0i,#2i #3  from  the  fourth  column  via 

9\  =  -arctan(j^) 

‘v 


02  =  -  arctanf  ■ 


#3  =  -arctanf' 


l2  +  l2 
*  T  y  N 

—  ' 

ni  +  ii  +  11 


Step  2:  Left-multiply  the  third  column  (kx,  ky,kz,kw)T  by  f?zw4(  -^3)-ffYZ4(-^2)^XY4(_^i  1 
to  transform  it  into  (kX3,ky3,kZ3,0)T ,  i.e.  the  body-fixed  vector  k  expressed  in  the 
third  coordinate  system.  Then  calculate  04,05  via 

9\  =  -  arctan( 

kv* 
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#5 


-  arctan( 


) 


Step  3:  Left-multiply  the  second  column  ( jx,jy,jz,jw)T  by 

#YZ3(  ~0s)#XY3(  -^4)-RzW4(  ~^3)-Ry24(  -^2  )^XY4(  -^1 ) 

to  transform  it  into  (;rs,  jys, 0,0)T,  i.e.  the  body-fixed  vector  j  expressed  in  the  fifth 
coordinate  system.  Then  calculate  via 

ffe  =  -  arctan(^-) 

Jys 

Note  that  steps  2  and  3  can  be  combined  into  one  by  3D  inverse  Euler  angle  formulas 
[17].  We  stress  again  the  similarity  between  the  rotation  chains  in  TZ4,TZ3,  TZ2.  Within  the 
rotation  chain  in  TZ'  (i  =  4,3,2),  only  the  angle  of  i?xYi  has  a  full  range  of  2x,  others 
are  restricted  in  a  range  of  x  to  eliminate  ambiguities.  That  is  why  there  is  no  i  before 
the  square  roots.  Figure  2  shows  the  body-fixed  vector  1  and  its  projections  in  TZ3  and  TZ2. 
The  angles  0j,02, d3  can  be  considered  as  polar  coordinates  extended  in  to  TZ4. 

Note  that  numerical  difficulties  arise  when  both  the  numerator  and  the  denominator  in 
the  argument  of  arctan  are  close  to  zero.  This  can  happen  when  at  least  one  of  the  angles 
02, 03,05  is  close  to  Jfcx  (k  =  0,  ±1,. . .),  and  is  a  well-known  drawback  of  Euler  angles  [17]. 

3.2  Projection  Centers  Specification 

We  specify  the  projection  centers  in  the  same  way  as  object  orientation.  Consider  the  screen 
for  display  as  a  rigid  body  and  attach  to  it  a  coordinate  system  (xe,y&,  z6,w6)  in  such  a 
way  that  the  first  projection  center,  called  the  sensor,  is  along  the  direction  of  w6-axis, 
and  the  second  projection  center,  called  the  eye,  is  along  the  direction  of  Z6-axis.  Using 
Euler  angles,  the  sensor's  direction  is  determined  by  01,02,03-  111  effect,  these  angles  are 
polar  coordinates  extended  in  TZ4.  The  eye’s  direction  is  determined  by  94,d 5,  which  are 
the  usual  polar  coordinates  in  3D,  also  called  the  azimuth  and  incidence  angles.  The  twist 
angle ,  06,  only  affects  the  orientation  of  the  final  2D  picture  in  the  projection  plane.  In 
addition  to  these  angles,  the  user  interface  offers  other  parameters  such  as  the  reciprocal  of 
the  distance  from  the  sensor  or  the  eye  to  the  origin,  the  field-of-view  angle,  the  distance 
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w 


Figure  2:  Body-fixed  vector  1  and  its  projections  in  TZ3  and  TZ? 

and  the  pitch  of  the  front  clipping  plane.  Ail  these  quantities  are  with  respect  to  the 
(xe,  ye,  Z6,  tee)-coordinate  system. 

4  Examples 

4.1  Hypersphere 

The  unit  hypersphere  S3  in  TZ4  has  the  equation: 

x2  +  y2  +  z2  +  w2  -  1  =  0 

If  the  hypersphere  is  projected  into  TZ3,  we  get  an  ellipsoid,  which  contains  little  information 
unless  there  is  a  4D  shading  model  resulting  in  different  intensities  at  different  interior  points 
of  the  ellipsoid.  A  more  convenient  way  to  visualize  the  hypersphere  is  to  mark  or  color 
certain  areas  on  it  in  analogy  to  drawing  longitude  and  latitude  circles  on  a  sphere  S2.  We 
can  obtain  the  “longitude”  and  “latitude”  on  a  hypersphere  by  the  Hopf  map  [11,  14]. 

h:S3  —  S2 
(x,y,z,w)  i — ♦  (X,Y,Z) 

where 

X  =  2(xy  +  zw) 


16 


Y  =  2 (xw  -  yz) 

Z  =  {x2  +  z2)  -  (y2  +  w2) 


For  example,  given  a  latitude  Z  =  a  on  52,  its  inverse  image  on  S3  is  a  2-torus,  S1  x  S1: 


x2  +  z 2 
y 2  +  w2 


1  +  a 
2 

1  -  a 
2 


0 

0 


In  Figures  6-8,  the  inverse  images  of  the  latitudes  Z  =  -  ^ ,  Z  =  0,  Z  =  ^  are  represented 
by  the  green,  gray  and  red  surfaces  respectively.  The  inverse  image  of  the  south  and  north 
poles,  Z  =  —  1  and  Z  =  1  are  represented  by  magenta  and  blue  curves.  We  also  display 
x,y,z,w  axes  as  straight  lines  of  red,  orange,  green,  and  blue,  respectively.  The  u;-axis  is 
now  invisible  because  the  sensor  is  positioned  on  it.  The  hypersphere  is  translated  by  1 
along  the  positive  to-axis  so  that  when  the  sensor’s  distance  is  2,  the  hypersphere  is  mapped 
to  the  whole  TZ3.  Figures  6-8  show  three  snap  shots  of  an  animation  where  the  sensor  is 
moving  from  infinity  towards  the  hypersphere.  The  moments  shown  are  when  the  sensor's 
distance  is  infinity,  3,  and  2,  respectively.  As  the  sensor  approaches,  the  circumference  of 
the  green  2-surface  enlarges  significantly  while  the  red  one  enlarges  only  moderately.  This 
animation  can  be  understood  by  comparing  it  with  an  animation  in  3D:  Look  at  Figure  6 
and  consider  the  surfaces  as  objects  in  3D  space.  Imagine  that  your  eye  is  moving  on  the 
y- axis  from  infinity  towards  the  origin,  then  you  will  see  a  similar  animation  but  happens 
in  3D  space.  Figure  8  also  shows  how,  by  using  the  tilted  front  clipping  plane,  some  hidden 
details  can  be  revealed.  See  also  [11]  for  a  visualization  of  4D  rotation. 


4.2  Offset  Curves 

Given  a  curve  f(x,y)  =  0  in  TZ2,  its  offset  curve  by  distance  r  >  0  can  be  formulated  by 
the  envelope  method  [5,  8]  as  a  set  of  equations: 

g  :  (i  -  u)2  +  (y  -  v)2  -  r2  =  0 

f(u,v)  =  0 

C:  Vuv9- 1  =  0 


where 


Vu  v9 


dg_  dg_  T 
[dudv) 
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t  -  MJ-l? 

~  {dv'  du ’ 


If  the  parametric  form  of  the  curve  /  is  available,  the  set  of  equations  can  be  simplified  as: 


h:  (x  -  u(t))2  +  (y  -  v(t))2  -  r2  =  0 

C' :  (x  -  u(t))u'(t)  +  (y  -  v(t))v'(t)  =  0 

Note  that  the  condition  C'  is  equivalent  to  =  0.  If  the  greatest  common  divisor  4>{t)  - 
GCD(u'(t),i/(<))  is  not  a  constant,  the  condition  C'  can  be  further  simplified  as  [5]: 


C"  :  (x  -  u(<))p(*)  +  (y  -  v(t))q(t)  =  0 


where 


p(t)  = 


u'(t) 


q(t)  = 


At) 


An  implicit  equation  for  the  offset  curve  can  be  determined  by  the  resultant  method  [5],  or 
using  Grobner  basis  [8].  The  offset  curve  can  also  be  traced  numerically  in  TZ4  or  7l3  by  the 
method  described  in  [2]. 

It  is  important  to  note  the  following  points  about  the  envelope  method  for  formulating 
offsets: 


1.  The  offset  curve  may  have  cusps  and/or  self-intersections  in  XY-plane  (see  Figure  3). 
But  the  singularities  often  disappear  when  the  curve  is  traced  in  higher  dimensional 
space. 

2.  The  equations  may  describe  additional  points  which  have  a  distance  r  from  the  singular 
points  on  the  curve  /  (see  Figure  4). 

We  tried  to  explain  these  phenomena  by  means  of  4D  surface  visualization.  The  equa¬ 
tions  <7  =  0,  /  =  0  are  two  hypersurfaces  in  XYUV-space  and  their  intersection  5  is  a 

2-surface.  Moreover,  at  the  point  p  =  (x,y,u,v)  on  S,  the  two  normals  are: 

nl  =  Vff(p)  =  (2(x  -  u),2(y  -  v),-2(x  -  u),-2(y  -  v))T 

n,  =  V/(p)  = 

They  are  linearly  independent  as  long  as  n2  is  a  nonzero  vector  since  (x  -  u)  and  (y  -  v) 
cannot  be  both  zero.  The  condition  C  can  be  rewritten  as  det(i,j,  nx,  n2)  =  0.  If  p 
is  a  nonsingular  point  on  S,  by  Lemma  3  (c)  it  is  a  silhouette  point  with  respect  to  an 
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orthographic  projection  with  two  centers  along  u-  and  v-axes.  The  silhouette  points  form 
a  curve  on  the  pipe-like  surface  5  in  U4.  In  Figure  9  and  10  we  show  the  2-surface  5 
and  the  silhouette  curve  corresponding  to  the  offset  curve  in  Figure  3.  In  1Z4  the  curve  is 
smooth  without  cusps  or  self-intersections  as  we  can  see  in  Figures  10  from  different  viewing 
direction. 

On  the  other  hand,  if  p  is  a  singular  point,  then  ni  and  n2  are  linearly  dependent,  and 
so  n2  must  be  a  zero  vector.  Surely  condition  C  is  satisfied,  but  according  to  our  definition 
they  are  not  silhouettes.  They  are  exactly  the  additional  points  described  in  the  second 
phenomenon.  In  Figure  11  and  12  the  2  surface  S  and  the  silhouette  curve  corresponding  to 
the  offset  curve  in  Figure  4  are  shown  in  different  viewing  directions.  The  silhouette  curves 
are  still  smooth  without  cusps  or  self-intersections.  But  the  2-surface  is  not  a  smooth  “pipe”. 
The  singular  points  form  a  circle  corresponding  to  the  dashed  circle  in  Figure  4. 

If  the  curve  /  has  a  parametric  form,  the  offset  curve  can  be  traced  in  XYT-space.  The 
two  equations  h  =  0,  =  0  are  two  surfaces  and  their  intersection  is  a  curve.  Note 

that  =  0  is  equivalent  to  y/i  •  k  =  0.  It  also  means  that  the  intersection  curve  is  the 
silhouette  on  the  surface  h  with  respect  to  a  orthographic  projection  along  the  t  axis.  But 
the  surface  h  =  0  is  smooth  without  any  singular  points  because  ||  and  cannot  be  zero 
simultaneously  and  so  sjh  is  always  a  nonzero  vector.  The  dashed  circle  in  Figure  2  is 
actually  another  branch  of  silhouette  curve  as  shown  in  Figures  13  and  14. 

If  the  greatest  common  divisor  <p(t)  is  not  a  constant,  the  condition  y/i  ■  k  =  0  is 
equivalent  to: 

<t>(t)[(x  -  u{t))p(t)  +  (y  -  v{t))q(t)\  =  0 

The  factor  4>(t)  —  0  represents  those  silhouette  curve  branches  that  are  circles  resulting 
from  intersecting  the  pipe-like  surface  h  =  0  with  the  planes  t  =  f,  perpendicular  to  the 
t-axis,  where  t,’s  are  the  roots  of  <p(t).  The  other  factor  is  the  same  as  condition  C",  and 
represents  the  silhouette  curve  branches  corresponding  to  the  offset  curve. 

4.3  Collision  Detection 

Cameron  [3]  proposes  an  algorithm  to  detect  collisions  based  on  4D  model.  The  basic  idea 
is  that  “if  an  object  can  be  represented  by  a  set-combination  in  a  CSG  scheme,  and  the 
primitive  objects  can  be  extruded  (into  4D)  in  this  scheme,  then  the  extrusion  of  the  object 
is  the  set-combination  of  the  extrusions  of  the  primitives”.  The  extruded  object  means  the 
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object  in  motion  considered  in  XYZV-space.  Two  moving  objects  collide  if  and  only  if  the 
intersection  of  their  extrusions  is  nonempty.  The  problem  is  then  reduced  to  testing  the 
intersection  of  each  pair  of  primitive  extrusions.  In  3D  the  primitives  can  be  the  half  space 
to  one  side  of  an  oriented  surface.  The  extrusion  of  the  surface  is  a  hypersurface  in  4D 
space. 

We  discuss  some  points  in  the  design  of  an  algorithm  to  test  whether  two  hypersurfaces 
intersect.  The  hypersurfaces  are 

f(x,y,z,t)  =  0,  g{x,  y,  z,  t)  =  0 

The  intersection  is  a  2-surface  in  T2.4  and  can  be  examined  with  our  visualization  system. 

Consider  a  cylinder  of  radius  rc  about  the  x-axis  moving  in  the  positive  y-direction 
at  a  constant  speed  vc,  and  a  sphere  of  radius  r3  moving  in  the  negative  x-direction  at 
a  constant  speed  vs.  At  the  time  t  =  0  both  are  at  the  origin  as  shown  by  the  dashed 
cylinder  and  sphere  in  Figure  5.  We  consider  the  2-surface  that  is  the  intersection  of  the 
two  hypersurfaces: 

(y  -  vct)2  +  z2  -  r]  =  0 
x2  +  y2  +  (z  +  v,t)2  -  r2  =  0 

Note  that  as  the  radii  and  speeds  change,  the  topology  of  the  resulting  intersection  2-surface 
can  be  quite  different.  Figure  15  shows  the  case  rc  =  r,  =  1,  vc  =  0,  v,  =  1.  The  green 
wire  mesh  represents  the  the  cylinder  at  rest,  i.e.  intersection  of  the  cylinder  with  the 
hyperplane  t  =  0.  Since  the  sensor's  position  is  just  a  little  off  the  t  axis,  the  intersection 
2-surface  resembles  the  sweep  of  the  intersection  curve  in  V? ■  Figure  16  shows  the  case 
rc  =  1,  rt  =  0.7,  vc  =  0.2,  v,  =  1.  Since  the  radius  of  the  sphere  is  smaller,  the  intersection 
2-surface  has  two  separate  components.  Figure  17  shows  the  case  rc  =  1,  r,  =  1.2,  vc  =  0. 
v,  =  1.  The  sphere  has  a  larger  radius.  Although  the  2-surface  is  connected,  it  has  a  “hole" 
due  to  the  fact  that  at  a  certain  time  period  the  sphere  and  the  cylinder  intersect  in  a  curve 
with  more  than  one  branch. 

One  advantage  of  using  4D  geometry  is  that  the  intersection  is  nonempty  as  long  as 
there  is  a  collision.  The  commonly  used  method  to  test  intersection  in  7 13,  at  a  series  of 
sampling  times  may  not  be  able  to  detect  the  following  cases: 

1.  The  interpenetration  begins  and  ends  between  two  consecutive  sampling  times. 
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V, 


r-0 


Figure  5:  A  cylinder  and  a  sphere  in  motion 

2.  If  in  TZ3  the  algorithm  does  only  boundary  intersection  test  but  not  containment  test, 
an  interpenetration  that  lasts  longer  than  one  sampling  period  may  also  be  overlooked. 

Often  we  only  need  to  find  the  initial  colliding  point ,  i.e.  the  point  in  the  intersection 
2-surface  with  the  smallest  value  of  t.  Since  we  do  not  consider  “patches”  or  “trimmed 
surface”,  a  necessary  condition  for  initial  colliding  point  is: 

df/dx  df/dy  df/dz 
dg/dx  dg/dy  dg/dz 

Together  with  /  =  0,  g  —  0  it  represents  a  zero-dimensional  point  set  in  II4.  If  /,  g  are 
polynomials,  this  set  can  be  solved  by  algebraic  methods,  e.g.  using  the  Grobner  Basis 
techniques  [8].  When  the  object  motion  includes  a  rotation,  /  and  g  are  no  longer  algebraic 
hypersurfaces.  Such  a  set  of  equations  often  is  hard  to  solve.  One  possible  method  is  to 
relax  the  condition  C  as: 

r, .  df_dg_  _  df_dg_  =  „ 

dx  dy  dy  dx 

Together  with  /  =  0,  g  =  0  it  represents  a  curve  in  TZ4  and  can  be  traced  numerically.  The 
point  on  the  curve  with  the  smallest  t  is  the  initial  colliding  point. 

The  condition  C'  can  be  rewritten  as: 

det(l,  k,  v/>  Vff)  =  0 
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According  to  Lemma  3  (c),  this  is  the  silhouette  curve  on  the  2-surface  with  respect  to  a 
orthographic  projection  along  x-  and  y-axes.  It  also  works  if  k  is  replaced  by  any  nonzero 
linear  combination  of  i,  j,  k,  corresponding  to  a  projection  in  other  directions.  It  is  important 
to  note  the  following  phenomena  when  tracing  the  silhouette  curve: 

1.  The  2-surface  may  have  several  separate  components.  Then  the  silhouette  curve  must 
have  several  separate  branches  as  seen  in  Figure  16. 

2.  If  the  2-surface  is  connected,  it  is  nevertheless  possible  that  the  silhouette  curve  has 
several  separate  branches.  See  also  Figure  17. 

Finding  the  initial  colliding  point  can  also  be  considered  as  a  constrained  nonlinear 
programming  problem  optimizing  t  (see,  e.g.  [12]).  Investigating  the  topology  of  the  2- 
surface  by  means  of  4D  visualization  may  help  to  design  faithful  algorithms. 

5  Conclusion 

We  have  discussed  two  aspects  of  visualization  of  2-surfaces  in  4D  space:  first,  some  geo¬ 
metric  issues  related  to  developing  an  interactive  4D  visualization  system,  and  second,  some 
examples  illustrating  the  use  of  such  a  system. 

We  have  presented  a  method  suitable  for  implementation  on  conventional  z-buffer  based 
graphics  workstations.  The  input  2-surfaces  are  polygonalized  in  4D  space.  The  4D  poly¬ 
gons,  with  each  vertex  attached  to  two  linearly  independent  normal  vectors,  are  then  pro¬ 
jected  into  3D  space  and  fed  into  the  3D  graphics  engine.  The  normal  of  the  projected 
surface  can  be  calculated  efficiently  by  Lemma  1.  The  main  problem  in  devising  a  user 
interface  is  how  to  specify  the  orientation  of  objects  and  the  directions  of  projection.  They 
can  be  done  uniformly  by  the  Euler  angles  when  extended  properly  into  the  4D  space. 

In  Lemma  2  and  3  we  have  also  discussed  the  conditions  for  a  point  on  the  2-surface  to  be 
a  silhouette  point  with  respect  to  a  projection.  This  concept  plays  an  important  role  in  the 
explanation  of  our  examples.  Offset  curves,  when  traced  in  4D  space,  can  be  considered  as 
silhouette  curves  of  a  pipe-like  2-surface.  Silhouette  curves  also  roughly  describe  the  shape 
of  the  2-surface  that  is  the  intersection  of  two  moving  surfaces  in  3D  space.  We  pointed 
out  some  phenomena  to  be  noticed  in  designing  faithful  algorithms  for  tracing  offset  curves 
and  searching  for  initial  colliding  points.  Our  experience  convinced  us  that  such  a  4D 
visualization  system  is  useful  for  understanding  four  dimensional  geometry. 
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Some  problems  remain  to  be  solved.  For  example,  the  quality  of  the  pictures  may  not 
be  good  enough  because  of  the  resolution  of  the  polygonalization.  When  the  sensor  or  the 
eye  is  too  close  to  the  2-surface,  this  problem  becomes  especially  acute,  since  the  nearby 
polygons  are  seen  much  larger  than  those  further  away. 
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Figure  6  Hypersphere  viewing  from  the  sensor  at 


2 


Figure  7  Hypersphere  viewing  from  the  sensor  at  distance  3 


Figure  9  Offset  curve  of  Figure  3  traced  in  4D,  viewing  orthographically 


of  Figure  3  traced  in  4D,  viewing  perspectively 


z-  and  w-axes  are  invisible 


Figure  1 1  Offset  curve  of  Figure  4  traced  in  4D,  viewing  orthographically 


y 


w-axis  is  hidden  by  the  2-surface 


Figure  12  Offset  curve  of  Figure  4  traced  in  4D,  viewing  perspectively 


z-axis  is  invisible 


Figure  13  Offset  curve  of  Figure  4  traced  in  3D,  viewing  orthographically 


Figure  14  Offset  curve  of  Figure  4  traced  in  3D,  viewing  perspectively 


Figure  15  Intersection  of  a  cylinder  and  a  moving  sphere  with  the  same  radius 


